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We report high statistics measurements of inclusive charged hadron production in Au+Au and 
p+p collisions at x /s NN =200 GeV. A large, approximately constant hadron suppression is observed 
in central Au+Au collisions for 5<pt<12 GeV/c. The collision energy dependence of the yields and 
the centrality and pr dependence of the suppression provide stringent constraints on theoretical 
models of suppression. Models incorporating initial-state gluon saturation or partonic energy loss 
in dense matter are largely consistent with observations. We observe no evidence of pr-dependent 
suppression, which may be expected from models incorporating jet attentuation in cold nuclear 
matter or scattering of fragmentation hadrons. 

PACS numbers: 25.75.Dw,25.75.-q,13.85.Hd 

High energy partons propagating through matter are charge density of the medium This process can pro- 
predicted to lose energy via induced gluon radiation, with vide a sensitive probe of the hot and dense matter gen- 
the total energy loss strongly dependent on the color erated early in ultrarelativistic nuclear collisions, when a 
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plasma of deconfined quarks and gluons may form. The 
hard scattering and subsequent fragmentation of partons 
generates jets of correlated hadrons. In nuclear collisions, 
jets may be studied via observables such as high trans- 
verse momentum (high px) hadronic inclusive spectra Q 
and correlations. Several striking high px phenomena 
have been observed at the Relativistic Heavy Ion Col- 
lider (RHIC) [ESS 00 including strong suppression 
of inclusive hadron production 0, 0,0 . These phenom- 
ena are consistent with large partonic energy loss in high 
energy density matter [H IslM Hoi Hl| though other mech- 
anisms have been proposed, including gluon saturation 
in the initial nuclear wavefunction 1 121 . attenuation of 
jet formation in cold nuclear matter |l3j . and scattering 
of fragmentation hadrons • Additional measurements 
are required to discriminate among these pictures and to 
isolate effects due to final state partonic energy loss. 

We report high statistics measurements by the STAR 
Collaboration of the inclusive charged hadron yield 
(h + + h~)/2 (defined as the summed yields of primary 
7r^~ , K ± , p and p) in Au+Au collisions and non-singly 
diffractive (NSD) p+p collisions at nucleon-nucleon cen- 
ter of mass energy v /s NN =200 GeV. The Au+Au data ex- 
tend considerably the pr range of earlier charged hadron 
suppression studies, and the p+p data are the first such 
measurement at this energy. Comparisons are made to 
several theoretical models. The high precision and broad 
kinematic coverage of the data significantly constrain the 
possible mechanisms of hadron suppression. In addition, 
the energy dependence of the yields may be sensitive to 
gluon shadowing at low Bjorken x in heavy nuclei. 

We compare the data to two calculations based on 
hard parton scattering evaluated via perturbative QCD 
(pQCD-I [3 and pQCD-II 9]) and to a calculation ex- 
tending the saturation model to high momentum trans- 
fer [l2j. Both pQCD models for Au+Au collisions in- 
corporate nuclear shadowing, the Cronin effect |lfij . and 
partonic energy loss, but with different formulations. 
pQCD-I results excluding one or more nuclear effects 
are also shown. Neither pQCD calculation includes 
non-perturbative effects that generate particle species- 
dependent differences for px<5 GeV/c|l.4 Il7j. 

Charged particle trajectories were measured in the 
Time Projection Chamber (TPC) 0. The magnetic 
field was 0.5 T, resulting in a factor of three improve- 
ment in pt resolution at high px relative to [!E1. After 
event selection cuts, the Au+Au dataset comprised 1.7 
million minimum bias events (97 ± 3% of the geomet- 
ric cross section cr^^ 1 ) and 1.5 million central events 
(10% of cr^om"). Centrality selection and analysis of 
spectra follow Ref. 0- Background at high px is dom- 
inated by weak decay products, with correction factors 
calculated using STAR measurements of A(+S°) and K® 
for pr<6 GeV/c and assuming constant yield ratios 
A(+E )/(h++h") and K°/(h++b.-) for p T >6 GeV/c. 
The A(+E°) yield was scaled by a factor 1.4 to account 



TABLE I: Multiplicative correction factors applied to the 
measured yields at pr=W GeV/c for p+p and Au+Au data. 
Factors vary by approximately 5% within 4<pr<12 GeV/c 
and have similar uncertainties. 





Tracking 


Background 


Pt resolution 


P+P 

Au+Au 60-80% 
Au+Au 0-5% 


1.18 ± 0.07 
1.11 ± 0.06 
1.25 ± 0.06 


0.90 ± 0.08 
0.95 ± 0.05 
0.94 ± 0.06 


0.89 ig;« 
0-97 

0-95 1°;° 5 5 



for £+ decays. Table[Q summarizes the correction factors 
at high px ■ 

After event selection cuts, the p+p dataset comprised 5 
million mainly NSD events, triggered on the coincidence 
of two Beam-Beam counters (BBCs). The BBCs are an- 
nular scintillator detectors situated ±3.5 m from the in- 
teraction region, covering pseudorapidity 3.3<| 77 1< 5.0. A 
van der Meer scan [5^ measured the BBC trigger cross 
section to be 26.1+0. 2(stat)±1.8(sys) mb. The BBC trig- 
ger was simulated using PYTHIA 21] and HERWIG 
[22j events passed though a GEANT detector model. 
The PYTHIA trigger cross section is 27 mb, consistent 
with measurement, of which 0.7 mb results from singly 
diffractive events. The PYTHIA and HERWIG simula- 
tions show that the trigger accepts 87+8% of all NSD 
events containing a TPC track, with negligible track px- 
dependence. Non-interaction backgrounds contributed 
3+2% of the trigger rate. The high p+p interaction 
rate generated significant pileup in the TPC. Valid tracks 
matched an in-time hit in the Central Trigger Barrel 
(CTB 18] ) surrounding the TPC and projected to a 
distance of closest approach DCA<\ cm to the aver- 
age beam trajectory. To avoid event selection multi- 
plicity bias, an approximate event vertex position along 
the beam (z ver t) was calculated by averaging zdca over 
all valid tracks. Accepted events were required to have 
I z V ert I < 75 cm, corresponding to 69+4% of all events. 
The track px fit did not include the event vertex. The 
CTB track-matching efficiency is 94+2% and combina- 
torial background is 2+2%. Other significant p+p track- 
ing backgrounds result from weak decays and antinucleon 
annihilation in detector material, with corrections calcu- 
lated using HIJING jll] and preliminary STAR measure- 
ments. Correction factors at high pt are given in Tabled 
For p+p collisions relative to peripheral Au+Au, exclu- 
sion of the event vertex from the pt fit results in poorer 
Pt resolution, while the CTB matching requirement re- 
sults in lower tracking efficiency. The p+p inclusive spec- 
trum was also analysed for px<3.5 GeV/c by an inde- 
pendent method in which a primary vertex is found and 
incorporated into the track fit, with consistent results. 

The p+p NSD differential cross section is the product 
of the measured per-event yield and the BBC NSD trig- 
ger cross section, and has a normalization uncertainty 
of ±14%. The charged hadron invariant cross section 
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FIG. 2: -R200/ 130 (Pt) vs. pr for (h + + h~)/2 for four differ- 
ent centrality bins. The overall normalization uncertainty is 
t.10% for the 40-60% bin and is negligible for the other panels. 
Calculations are described in the text. 



FIG. 1: Inclusive invariant pr distributions of (h + + h~ )/2 for 
centrality-selected Au+Au and p+p NSD interactions. Hash 
marks at the top indicate bin boundaries for pt>4 GeV/c. 
The invariant cross section for p+p is indicated on the right 
vertical axis. 



has been measured in p + p collisions at y / s=200 GeV 
[24l | . The p+p cross section reported here is smaller by a 
factor of 0.79 ± 0.18, approximately independent of pr, 
where the uncertainty includes the two spectrum normal- 
izations and the correction for different acceptances [3j. 
The difference is due in large part to differing NSD cross 
section, which is 35 ± 1 mb in [24| but is measured here 
to be 30.0 + 3.5 mb. 

Figure ^ shows inclusive invariant px distributions of 
(h+ + h~)/2 within |?/|<0.5 for Au+Au and p+p colli- 
sions at ^s NN =200 GeV. The Au+Au spectra are shown 
for percentiles of o^om 1 with 0-5% indicating the most 
central (head-on) collisions. Error bars are the quadra- 
ture sum of the statistical and systematic uncertainties 
and are dominated by the latter except at the highest px- 

Figure shows -R20o/13o(pt), the ratio of charged 
hadron yields at ^/s NN =200 and 130 GeV0], for central- 
ity selected Au+Au collisions. Error bars are the quadra- 
ture sum of the statistical and systematic uncertainties, 
dominated for pr>4 GeV/c by statistics at 130 GeV. 
In the absence of nuclear effects, the hard process in- 
clusive yield in nuclear collisions is expected to scale as 
(Nbin), the average number of binary collisions for the re- 
spective centrality selection. R2oo/13o(pt) has not been 



scaled by the ratio N h; n ( 200) /Nhin (130), which Glauber 
model calculations pi l2fij give as ~ 1.02 for all centrali- 
zes. Figurc|21also shows the saturation model calculation 
and pQCD-I calculations for p+p and centrality-selected 
Au+Au collisions (shadowing-only and full). Both mod- 



els approximately reproduce the pr-dependence of the 
ratio for Au+Au for p T >2 GeV/c, with pQCD-I slightly 
better for more peripheral collisions. The various pQCD- 
I calculations shown illustrate that in this model the 
reduction in i?20o/i3o(PT) f° r Au+Au relative to p+p 
is predominantly due to nuclear shadowing|l5j. This 
sensitivity arises because the shadowing is a;-dependent 
and at fixed px, different y/s corresponds to different 
xt = 2pt/\/s. The quantitative agreement of pQCD- 
I with the data improves for more peripheral collisions, 
suggesting that the prescription for the centrality depen- 
dence of shadowing in |l5j may not be optimal. Alterna- 
tively, introduction of y^-dependent energy loss to the 
model in |15| may also improve the agreement. 

Nuclear effects on the inclusive spectra are measured 
by comparison to a nucleon-nucleon (NN) reference via 
the nuclear modification factor: 



Raa(pt) 



d 2 N AA /dp T dr] 
T AA d 2 a NN /d PT df] 



(1) 



where T , A J 4=(Nbm)/ <7 ^^ from a Glauber calculation ac- 
counts for the nuclear collision geometry 0, [2{| and we 
adopt cr^ei = 42 mb. d 2 <r NN j dp T dr\ refers to inelas- 
tic collisions, whereas we have measured the p+p NSD 
differential cross section. However, singly diffractive in- 
teractions contribute predominantly to low px H3- A 
multiplicative correction based on PYTHIA, applied to 
d 2 a NN /dp T dr] in Eq. □ is 1.05 at p T =0.4 GeV/c and 
unity above 1.2 GeV/c. 

Figure |3] shows Raa(pt) at ^/s NN =200 GeV for 
centrality-selected Au+Au spectra relative to the mea- 
sured p+p spectrum. Horizontal dashed lines show 
Glauber model expectations 0, H^J for scaling of the 
yield with (Nbin) or mean number of participants (N part ), 
with grey bands showing their respective uncertainties 
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FIG. 3: Raa(pt) (Eq. □ for (h + +h _ )/2 in |r/|<0.5, for 
centrality-selected Au+Au spectra relative to the measured 
p+p spectrum. The p+p spectrum is common to all panels. 
Calculations are described in the text. 



summed in quadrature with the p+p normalization un- 
certainty. The error bars represent the quadrature sum of 
the Au+Au and remaining p+p spectrum uncertainties. 
For pt<6 GeV/c, Raa{pt) is similar to that observed 
at ^/s NN =130 GeV [3j. Hadron production for 6<px<10 
GeV/c is suppressed by a factor of 4-5 in central Au+Au 
relative to p+p collisions. 

Figure|3]also shows the full pQCD-I calculation and the 
influence of each nuclear effect. The energy loss for cen- 
tral collisions is a fit parameter, with the px and central- 
ity dependence of the suppression constrained by theory. 
The Cronin enhancement and shadowing alone cannot 
account for the suppression, which is reproduced only if 
partonic energy loss in dense matter is included. The full 
calculation generally agrees with data for px>5 GeV/c if 
the initial parton density in central collisions is adjusted 
to be ~15 times that of cold nuclear matter j^. pQCD- 
II exhibits similar agreement for central collisions. In 
Ref. 0, the pQCD-II calculation was used to predict a 
PT-independent suppression factor in this px range from 
the interplay between shadowing, the Cronin effect, and 
partonic energy loss. 

Figure 0] shows Rcp(pt), the (Nbin)-normalized ratio 
of central and peripheral Au+Au spectra. Rcp{pt) ex- 
tends to higher px than Raa(pt), with smaller point- 
to-point uncertainties. The error bars show the quadra- 
ture sum of statistical and uncorrelated systematic un- 
certainties. Statistical errors dominate the uncertainties 
for pt>8 GeV/c. Lines and grey bands are as in Fig. [3] 
Rcp(pt) for pt<6 GeV/c is similar to measurements 



0.1 



o 
or. 



0.1 



Binary Collision scaling 



Participant scaling 
0-5%/40-60% 



pQCD-l 
pQCD-II 
Saturation 



0-5%/60-80% 



FIG. 4: R C p(pt) vs. p T for (h 4 
described in the text. 



6 8 10 12 
p T (GeV/c) 

h~)/2. Calculations are 



at ^/s NN =130 GeV 3], but is now seen to be approxi- 
mately constant for 5<pt<12 GeV/c. It is consistent 
with (Npart) scaling at px ~ 4 GeV/c as reported in 
but is below (N par t) scaling at higher px- 

The pT-dependence of the suppression in Figure 0] is 
well reproduced for px>5 GeV/c by the full pQCD-I and 
pQCD-II calculations in both panels and the saturation 
calculation in the upper but not the lower panel. The 
magnitude of suppression is fitted to the central colli- 
sion data in the pQCD models but is predicted in the 
saturation calculation. Attenuation of initial jet forma- 
tion due to multiple nucleon interactions [13j generates 
an increase in partonic Raa{pt) for central collisions of 
a factor ~ 2 in h<Er<12 GeV. A similar pT-dependencc 
would be expected for high px hadrons, in contrast to 
observations. Suppression in the final state due to in- 
medium scattering of fragmentation hadrons should also 
result in a rising Raa(pt) with increasing pt due to the 
dependence of hadron formation time on the total jet en- 
ergy though detailed comparison of this model to 
data requires further theoretical development. 

In summary, STAR has measured inclusive charged 
hadron yields from Au+Au and p+p collisions at 
^/s NN =200 GeV, at higher precision and over a much 
broader px range than previous measurements. Large, 
constant hadron suppression is observed in central nu- 
clear collisions at high px- The systematic behaviour 
of the suppression at high px is well described both by 
pQCD calculations incorporating final-state partonic en- 
ergy loss in dense matter and a model of initial-state 
gluon saturation, though the latter model provides a 
poorer description of peripheral collision data. The isola- 
tion of initial state effects on high px hadron production 
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may be achieved through the study of d+Au collisions at 
RHIC, allowing a quantitative measurement of final state 
effects from the data presented here. 
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